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Abstract. A new model atmosphere analysis for the peculiar DQ white dwarf discovered by Carollo et al. (2002) is 
presented. The effective temperature and carbon abundance have been estimated by fitting both the photometric 
data (UBjVRfInJHK) and a low resolution spectrum (3500 < A < 7500 A) with a new model grid for helium-rich 
white dwarfs with traces of carbon (DQ stars). We estimate T c g ~ 5120 ±200 K and log[C/He]~ —5.8 ±0.5, which 
make GSC2U J131147.2±292348 the coolest DQ star ever observed. This result indicates that the hypothetical 
transition from C2 to C2H molecules around T e g = 6000 K, which was inferred to explain the absence of DQ stars 
at lower temperatures, needs to be reconsidered. 
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1. Introduction 

The presence of carbon in the atmospheres of some non- 
DA white dwarfs (defined as spectral type DQ) is gener- 
ally explained by the convective dredge-up from the stel- 
lar core to the outer photospheric layers (Koester et al. 
1982; Pelletier et al. 1986). C2 molecules are responsi- 
' ble for the absorption bands (e.g. in particular the Swan 
bands) which are the typical signature of the DQ stars. 
The spectral energy distribution of these stars changes 
significantly as a function of the effective temperature, 
T e ff, and carbon abundance, [C/He], as shown by the the- 
oretical atmosphere models of Koester et al. (1982) and 
Wegner & Yackovich (1984). Typically, strong absorption 
bands are expected for the coolest DQ stars, even in the 
case of low carbon abundances. 

However, past surveys revealed DQ stars with effective 
temperature above 6500 K only (Bergeron et al. 1997). 
The existence of this cut-off is not well understood. In fact, 
if cool DQ stars with T e g < 6500 K do exist, their strong 
Swan bands would result in peculiar colors and spectra 
which should make these objects easily recognizable. 

On the other hand, at low temperature carbon can be 
present also in a different form, as C2H molecules, if some 
hydrogen is also present. The electronic transition spectra 
of the C2H are not known from theory or laboratory ex- 
periments, but the observed spectra of the few known C2H 
stars show molecular absorption bands similar to the Swan 



bands shifted by about 150 A toward the blue. This shift 
cannot be explained as an effect of pressure shift of the 
Swan bands in a helium dominated atmosphere (Bergeron 
at al. 1994) or as a displacement due to a magnetic field 
(Schmidt et al. 1995). The presence of a certain fraction of 
hydrogen in the atmosphere of such non-DA white dwarfs 
can also be inferred by the collision induced absorption in 
the near IR due to H2 molecules, as in the case of LHS 
1126. 

These observations suggest the hypothesis that DQ 
white dwarfs turn into C2H stars when T e g is below 6500 
K, due to a not well identified physical mechanism that 
should inject hydrogen 1 in the He-dominated atmosphere 
of the DQ stars (Bergeron et al. 1997, 2001). 

Recently, Carollo et al. (2002) discovered 
GSC2UJ131147.2±292348 during a proper motion 
survey for halo white dwarfs based on the photographic 
material used for the construction of the GSC-11 (McLean 
et al. 2000). As shown in Figure 2 of Carollo et al. (2002), 

1 The evolution of the "missing" cool DQ stars is related to 
the more general problem of the non-DA gap which derives 
horn the apparent lack of non-DA stars observed with tem- 
perature 5100 ^ T c ff ^ 6100 K. At the moment, the cause of 
this effect, which could depend on the physical and chemical 
evolution of the white dwarf atmospheres during the cooling 
phases as well as on not sufficiently understood input physics, 
is not well established (see e.g. Bergeron et al. 1997, Malo et 
al. 1999). 
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Fig. 1. Two-color diagram constructed from VJK. For 
the theoretical predictions the continuous lines arc lines 
of constant carbon abundance from log[C/He] = —8 (top) 
to —4 (bottom); dotted lines are lines of constant effective 
temperature as described in the text. The cross marks 
GSC2U J131147.2+292348 which is the coolest object near 
5100 K, while the circles indicate other DQ stars. 

this object appears as a very peculiar carbon rich white 
dwarf due to the simultaneous presence of strong C2 
Deslandres-d'Azambuja and Swan bands, with an evident 
depression of the continuum in the Swan region between 
4500 and 6200 A. No other DQ star shows both these 
extreme features. 

GSC2UJ131147.2+292348 represents an enigma as 
well as an opportunity to test the predictions of the sce- 
nario discussed above for these objects. In fact, tempera- 
tures as low as 6000 K were already estimated by Carollo 
et al. (2002) by means of a simple spectral analysis. As de- 
scribed in the next sections, here we adopt a more sophis- 
ticated fitting technique based on new atmosphere models 
in order to estimate accurately the temperature and com- 
position of this object. The new results confirm cool values 
of the temperature, T c g ~ 5100 K, and the implications 
of this result will be briefly discussed. 

2. Model atmosphere analysis of 
GSC2U J131147.2+292348 

GSC2U J131147.2+292348 shows extremely strong bands 
of the C2 molecule, especially the Swan and Deslandres- 
d'Azambuja band systems in the optical range. As these 
bands obviously block a significant fraction of the total 
flux, they will influence the temperature structure of the 
atmosphere models. We have therefore calculated a new 
grid of model atmospheres, which takes into account the 
blanketing effect of 537 bands from the Swan, Deslandres- 
d'Azambuja, Fox-Herzberg, Phillips, and Ballik-Ramsay 



systems, as well as numerous atomic carbon lines and the 
resonance lines of Hel in the EUV. The general molecu- 
lar data were taken from Hubcr & Herzberg (1979), the 
Franck-Condon factors for the vibrational transitions from 
various sources (Krishnaswamy & Odell 1977; Dwivedi et 
al. 1978; Spindler 1965; Sharp 1983). Molecular absorption 
is treated with the "just overlapping line approximation 
(JOLA)" in the version as described in Zeidler-K.T. & 
Koester (1982). 

The general procedures and input physics of the model 
atmosphere calculations are very similar to the description 
in Finley et al. (1997). As it is practically impossible for 
very cool white dwarfs to determine effective temperature, 
surface gravity, and in our case the carbon abundance si- 
multaneously, we have held logg fixed at the canonical 
value of 8.00. T cff for the grid ranges from 10000 to 4600 K, 
the abundance ratio log[C/He] by numbers from —8 to —4. 

2.1. Magnitudes and colors 

For cool white dwarfs magnitudes and colors, especially 
in the infrared, are very useful for the determination of 
atmospheric parameters, a method pioneered by Bergeron 
et al. (1997). 

Since magnitudes in both the standard and pho- 
tographic system, with a spectral coverage from the 
ultraviolet to the near IR, have been observed for 
GSC2U J131147. 2+292348 , we have also calculated the- 
oretical magnitudes for our model grid. We adopted U,V 
from the photographic photometry given by Moreau & 
Reboul (1995), while Carollo et al. (2002) provided pho- 
tographic Bj, Hp and I^v in the natural photographic sys- 
tem of the POSS-II plates, plus standard JHK photome- 
try from observations carried out at the 4-m TNG (La 
Palma). The methods and the magnitude zeropoints used 
for the UV-JHK bandpasses in the standard Johnson sys- 
tem are described in detail in Zuckerman et al. (2003). 
For the photographic Bj, Hp and Ijv (approximately cor- 
responding to the Johnson-Cousins B(RI) C ) we have used 
the same transmission curves adopted for the photometric 
calibration of the GSC-II plates and determined the zero- 
points from integrations over the Vega flux as obtained 
from the STScI archive. 

The available magnitudes from U to K com- 
pletely determine the energy distribution of 
GSC2U J131147. 2+292348 . We have used our automatic 
least squares fitting routine, described in Zuckerman 
et al. (2003) to determine the best fitting parameters 
within the T c s - log[C/He] grid, resulting in an extremely 
low effective temperature around 5000 K (Table 0). The 
first row in the Table gives the observed magnitudes, 
the second the assumed errors. The third row are the 
theoretical predictions for the best fit parameters T e g- = 
4980 K, log[C/He] = -6.17. 

As an internal check, we tested the effect of fitting 
the physical parameters with only the standard UVJHK 
photometry with respect to the global solution including 
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Table 1. Observed standard and photographic magnitudes for GSC2U J131147. 2+292348 and theoretical fits. 
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Fig. 2. Observed spectrum of GSC2U J131147. 2+292348 (thin line) and theoretical model (thick line). 



also the Bj, Rf and Ijv magnitudes. The last two rows of 
Tabled show that the parameters change only very little, 
indicating that the GSC-II magnitudes are certainly very 
consistent with the overall energy distribution. 

Figure ^ shows the position of 
GSC2UJ131147.2+292348 in a special two-color di- 
agram J — V vs. J — K, using only standard magnitudes 
to be able to compare with other known DQ white 
dwarfs. Continuous lines are lines of constant carbon 
abundance, from -4.0 to -8.0, dotted lines are lines of 
constant effective temperature from 4600 K to 10000 K 
in steps of 200 K. As can be seen, this diagram is not 
very useful at temperatures above 7000 K, because of 
the competing direct effect of flux blocking and the 
indirect blanketing effect on the temperature structure. 
However, in the range 4600 - 6600 K and the abundances 
considered here, the diagram gives a clear indication 
of the atmospheric parameters. The cross at the lowest 



temperatures is GSC2U J131147.2+292348 , for which 
we would determine T e ff= 5100, log[C/Hc] = —6.0 from 
this position. The other 11 circles are observations of DQ 
white dwarfs from Bergeron et al. (1997) and Bergeron et 
al. (2001), which clearly are all much hotter, in agreement 
with temperatures derived in Bergeron et al. (1997). 

2.2. Spectral fitting 

The spectrum of GSC2U J131147. 2+292348 has been 
described in detail in Carollo et al. (2002). They con- 
cluded that the extremely strong bands of the Swan and 
Deslandres-d'Azambuja systems in the optical range are 
compatible with models calculated in Wegner & Yackovich 
(1984), whereas the energy distribution in the infrared 
could be explained by a blackbody distribution of around 
6000 K. With our consistent model atmospheres available, 
we can apply our standard spectral fitting technique (e.g. 
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Kocstcr ct al. 2001) with a Levenberg-Marquard algorithm 
(Press et al. 1992) to find the minimum \ 2 solution, using 
T e ff and log[C/He] as two free fitting parameters instead 
of the usual T e g and log g in the case of DA or DB white 
dwarfs. The quasi continuum was forced to fit the model at 
two positions (around 4150 and 7000 A), allowing for re- 
maining small calibration errors of the spectral flux. The 
resulting parameters for the best fit are T c g— 5200 K, 
log[C/He] = —5.53. Figure 2 shows the observed spec- 
trum together with the theoretical model corresponding 
to these parameters. Qualitatively, the theoretical model 
describes the main features of the spectrum, in particular 
the very strong band systems. In the details discrepancies 
remain, which may have a number of origins: the tempera- 
ture structure of the models, the equation of state in these 
very high pressure atmospheres, and, most likely, miss- 
ing bands, due to unkown Franck-Condon factors for the 
bands with highly excited lower levels, which are weak at 
laboratory conditions, but may be important in the much 
hotter stellar atmosphere. Nevertheless, we consider the 
fit satisfactory and a confirmation of the low temperatures 
derived from the photometry. 

2.3. Results and conclusion 

The fitting procedure for the photometry as well as for 
the spectrum provides formal errors, derived from the as- 
sumed statistical errors of the observations. These are very 
small — typically 30-40 K for T off and 0.05 for log[C/He] 
— and definitely unrealistic, because the errors are dom- 
inated by systematic errors of the models and reductions. 
These errors can be estimated only very roughly, taking 
the differences between the solutions from photometry and 
spectrum as a guide. Since we believe that the spectral re- 
sult is more reliable, we give it double weight and take 
as the final result for the atmospheric parameters T c g = 
5120±200 K and log[C/He] = -5.8±0.5. The distance mod- 
ulus obtained from the photometric solution is 3.69 mag, 
corresponding to a distance of 55 pc and to a tangen- 
tial velocity Kan = 4.74 • \xd ~ 125 km s" 1 (p = 0.48" 
yr -1 ). Adopting the same kinematics assumptions as in 
Carollo et al. (2002), we obtain galactic velocity com- 
ponents with respect to the LRS, (U,V) ~ (—115,-1) 
km s _1 . These values are well consistent with the velocity 
ellipsoid of the galactic halo (la) and are still consistent 
with the thick disk kinematics (2a), while the member- 
ship of GSC2U J131147.2+292348 to the thin disk appears 
much less probable. 

However, one needs to keep in mind that the results 
have been derived using a fixed surface gravity of log g= 
8.00. While we do not expect the atmospheric parameters 
to change much with log g, the distance modulus depends 
of course on the radius of the star, which depends strongly 
on the assumed surface gravity. Allowing for a plausible 
range of 7.5 - 8.5, the radius could be different up to ±30%, 
with the same change resulting for the distance and ve- 
locity. 



With T e ff about 5100 K this star is by far the coolest 
known "normal" DQ object. It is below the cutoff seen 
by Bergeron et al. (1997) near 6500 K and also below or 
at least at the lower edge of the so-called non-DA gap. 
It cannot be true therefore that all DQ turn into C2H 
stars when they cool down, and one obvious explanation 
could be that some stars completely avoid any accretion of 
hydrogen, which is a prerequisite for the formation of this 
molecule. However, the final explanation of this puzzle as 
well as others concerning the non-DA gap will likely need 
to wait for the discovery of more similar objects from the 
ongoing large scale survey like the SDSS. 
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